Black Hole-Neutron Star Mergers: Disk Mass Predictions 
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Determining the final result of black hole-neutron star mergers, and in particular the amount of matter remain- 
ing outside the black hole at late times, has been one of the main motivations behind the numerical simulation 
of these systems. Black hole-neutron star binaries are amongst the most likely progenitors of short gamma-ray 
bursts — as long as they result in the formation of massive (at least ~ 0.1M©) accretion disks around the 
black hole. Whether this actually happens strongly depends on the physical characteristics of the system, and 
in particular on the mass ratio, the spin of the black hole, and the radius of the neutron star. We present here 
a simple two-parameter model, fitted to existing numerical results, for the determination of the mass remain- 
ing outside the black hole a few milliseconds after a black hole-neutron star merger. This model predicts the 
remnant mass within a few percents of the mass of the neutron star, at least for remnant masses up to 20% of 
the neutron star mass. Results across the range of parameters deemed to be the most likely astrophysically are 
presented here. We find that, for IOMq black holes, massive disks are only possible for fairly large neutron 
stars (R NS > 12km), or quasi-extremal black hole spins (obh/Mbh > 0.9). We also use our model to discuss 
how the equation of state of the neutron star affects the final remnant, and the strong influence that this can have 
on the rate of short gamma-ray bursts produced by black hole-neutron star mergers. 

PACS numbers: 04.25.dg, 98.70.Rz, 04.40.Dg 
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I. INTRODUCTION 

The potential of black hole-neutron star (BHNS) mergers 
as progenitors of short gamma-ray bursts (SGRBs) and their 
importance as sources of gravitational waves detectable by 
ground-based interferometers such as Advanced LIGO and 
VIRGO d-d] have driven most recent studies of these sys- 
tems. Gamma-ray bursts, in particular, are a likely result if 
the neutron star is tidally disrupted, and the final outcome of 
the merger is a massive accretion disk around the black hole 
(see 0] and references therein). 

Numerical simulations have taught us that BHNS merg- 
ers can be divided into two broad categories: either tidal ef- 
fects are strong enough that the neutron star is disrupted be- 
fore reaching the innermost stable circular orbit (ISCO) of the 
black hole, or the neutron star plunges into the hole before 
tidal disruption occurs. In the first case, some material from 
the disrupted star remains outside the black hole for long peri- 
ods of time (~ 0. 1 — Is) in the form of an accretion disk, a tidal 
tail, and/or unbound ejecta. In the second case, however, the 
entire neutron star is rapidly accreted onto the black hole. To 
first order, the most important parameters determining the out- 
come of a BHNS merger are the mass ratio of the binary JD- 7 ], 
the spin magnitude of the black hole Ja-lal, its orientation H, 
and the size of the neutron star JH |9[ 1 1 CXI - The formation of 
massive accretion disks is more likely to occur for black holes 
of low mass (at least down to mass ratios A/bh/A/ns ~ 3) 
and high spins, and for large neutron stars. 

Studying these mergers is a complex problem, and accu- 
rate results can only be obtained through numerical simula- 
tions in a general relativistic framework: results using approx- 
imate treatments of gravity can lead to qualitative differences 
in the dynamics of the merger, and large errors in the mass of 
the final accretion disk or of any unbound material. Unfor- 
tunately, general relativistic simulations are computationally 
expensive, and only ~ 50 BHNS mergers have been studied 
so far (see JUIH for reviews of these results). Additionally, 



a majority of these simulations considered binaries with mass 
ratios A/bh/A/ns ~ 2 — 3, while population synthesis models 
indicate that mass ratios A/bh/A/ns > 5 are astrophysically 
more likely Il3l,[l4tl . Existing general relativistic simulations 
are also fairly limited in the physical effects considered: only 
a few include magnetic fields IU5ill6ll or nuclear theory -based 
equations of state ifioll . and none have considered neutrino 
emission (although neutrinos have been included in simula- 
tions of neutron star-neutron star mergers lfl7ll ). Magnetic 
fields and neutrinos do not affect significantly the disruption 
of the star, or the amount of matter remaining outside the 
black hole after merger, but they are critical to the evolution 
of the post-merger remnant, and to the modelling of electro- 
magnetic and neutrino counterparts to the gravitational wave 
signal emitted by black hole-neutron star mergers. 

Given the size of the parameter space to explore and the 
cost of numerical simulations, obtaining accurate predictions 
for the final state of the system for all possible configurations 
is only feasible through the construction of a model which ef- 
fectively interpolates between known numerical results. Such 
a model can also be of great help to determine which binary 
parameters should be used in numerical simulations in order 
to study a specific physical effect (e.g. massive disks) with- 
out having to run many different configurations. Taniguchi et 
al. [18] obtained a first estimate of the limit between configu- 
rations in which no disk is formed and those for which some 
matter remains at late time by studying quasi-equilibrium con- 
figurations of non-spinning binaries. More recently, Pannarale 
et al. Il9ll computed estimates for the mass remaining outside 
the black hole at late times through the use of a toy-model 
studying the tidal forces acting on the neutron star, repre- 
sented by a tri-axial ellipsoid. However, their model was fitted 
to numerical simulations which underestimated the remnant 
masses, and at a time when general relativistic results were 
only available at low mass ratios A/bh/A/ns ~ 2 — 3. The 
qualitative dependence of the remnant mass in the parameters 
of the binary is captured by their model, but the quantitative 
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results do not match more recent simulations, particularly at 
higher mass ratios J3l • 

In this paper we show that much simpler models comparing 
the estimated separation at which tidal disruption of the neu- 
tron star occurs (dtidai) and the radius of the ISCO (i?isco) 
can accurately predict the mass remaining outside the black 
hole at late times. We fit two such models (with different ap- 
proximations for dtidai) to a set of 26 recent numerical simu- 
lations covering mass ratios in the range A/bh/A/ns = 3 — 7, 
black hole spins up to 0.9 and neutron star radii 7?ns ~ 
11 — 16 km. The case of black hole spins misaligned with 
the orbital angular momentum is not considered here, and we 
limit ourselves to low eccentricity orbits (high eccentricities 
only occur when the binary is formed through dynamical cap- 
ture, e.g. in nuclear or globular clusters). Both models match 
the simulation results within their expected numerical errors, 
a few percents of the original mass of the neutron star. In 
Sec. in we describe the models used, and their physical in- 
spiration. Sec. [Ill] summarizes the numerical results used to 
calibrate the models, while Sec.|lV]gives the best-fit parame- 
ters, and discuss the quality of the fits. Finally, in Sec.[V] we 
show predictions of these models across the entire parameter 
space. We also discuss their strong dependence in the size of 
the neutron star, and potential implications for the rate of short 
gamma-ray bursts originating from BHNS mergers. 



II. TIDAL DISRUPTION MODELS 

The models used here to estimate the mass remaining out- 
side the black hole at late times are based on a compari- 
son between the binary separation at which tidal forces be- 
come strong enough to disrupt the star, dtidai, and the radius 
of the innermost stable circular orbit -Risco- Intuitively, if 
dtidai ?5 -Risco, the neutron star will plunge directly into 
the black hole and no mass will remain outside the hole af- 
ter merger. On the other hand, if dtidai ^ -Risco, the star will 
be disrupted. Some disrupted material will then form an ac- 
cretion disk, some will be ejected in a tidal tail and fall back 
on the disk over timescales long with respect to the duration 
of the merger (which typically lasts a few milliseconds), and 
some might be unbound. 

The separation dtidai at which tidal disruption occurs can be 
estimated in Newtonian theory by balancing the gravitational 
acceleration due to the star with the tidal accelaration due to 
the black hole: 
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where _Rns is the radius of the neutron star, A/ns an d A/bh 
are the masses of the compact objects, and we work in units 
in which G = c = 1. In general relativity, these quantities are 
not uniquely defined. In practice we will use the radius of the 
star in Schwarzschild coordinates and the ADM mass of the 
compact objects, all measured at infinite separation. 



As for the radius of the ISCO, it is given by |2j 

Zi = 1 + (1-xIh) 1/3 [( 1 + Xbh) 1/3 + (1-Xbh) 1/3 
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= 3 + Z 2 - sign(xBH) ^(3 - Zi)(3 + Z 1 + 2Z 2 ) (3) 



where xbh = obh/A/bh is the dimensionless spin parameter 
of the black hole. 

To construct a model for the fraction of the baryon mass 
of the star remaining outside the black hole at late times, we 
assume that this mass is entirely determined by the relative lo- 
cation of i?isco an d dtidai, in units of the neutron star radius. 
A first guess for the remnant mass M^ odel is then the linear 
model: 
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where a , (3° and 70 are the free parameters of the model, 
and Af^ s is the baryon mass of the neutron star. However, 
this simple prescription fits the numerical data rather poorly. 
In particular, it stronly underestimates the impact of the neu- 
tron star compactness Cns = Afpjs / -^ns on the result. This 
problem is not overly surprising: dtidai was derived in New- 
tonian gravity, but applied to compact objects. In particular, 
it predicts a finite radius for tidal disruption even if we re- 
place the neutron star by a non-spinning black hole (for which 
C = 0.5). To improve the model, we use instead a corrected 
estimate of the distance for tidal disruption, in which compact 
objects are more strongly bound: 

dtidai = dtidai (1 — 2Cns)- (5) 

This leads to the following model for the mass remaining out- 
side the black hole at late times, M 
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We could have added a constant term 7 as in 



Eq.|4] but find that this does not improve the quality of the fit. 
At the current level of accuracy of numerical simulations, we 
will show in Sec. [IV] that this simple model is in agreement 
with known results. 

A potential improvement on the model described by Eq. [6] 
is to compute the tidal effects from the Kerr metric instead 
of the Newtonian formula. Fishbone 1I2H1 obtained analytical 
results for these effects. Using his results leads to a correction 
to the value of the separation at which tidal disruption occurs: 
Ctidal = dtidai / Rns is then solution of the implicit equation 
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with k = A1bh/^?ns- We can then write the corrected model 
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In practice, M™^ cl gives results consistent with the simpler 



model M- de , 
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III. NUMERICAL RESULTS 

To fit the parameters a and /3 of our model, we consider 
recent results from numerical relativity in the range q = 3 — 
7, Xbh = - 0.9 and C NS = 0.13 - 0.18. We neglect 
simulations at lower mass ratios, which are astrophysically 
less likely and cannot be modelled accurately by the simple 
formula assumed here. Larger spins and more compact stars 
would be interesting to consider, but have yet to be simulated. 
Indeed, according to Hebeler et al. ll22ll . neutron stars of mass 
Mns ~ 1.4AI© could be in the range Cns = 0.15 - 0.22. For 
the same neutron star mass, Steiner et al. ll23"ll found the most 
likely compactness to be Cns = 0.17 — 0.19. 

We also limit the model to spins aligned with the orbital an- 
gular momentum and to low-eccentricity orbits. Misaligned 
spins have only been studied for one set of binary parame- 
ters 1 8], so that we do not have enough information about their 
influence on the disk mass to include them in the model. High- 
eccentricity mergers have been studied by East et al. I24II25I1 . 
As for misaligned black hole spins, however, the data does 
not cover enough of the parameter space to be included in 
our fit. Additionally, eccentricity is only an issue for binaries 
formed in clusters: field binaries are expected to have negli- 
gible eccentricities at the time of merger. Finally, we neg lect 
the influence of magnetic fields, as both Etienne et al. IU5I1 and 
Chawla et al. |[l6ll find their effect on the remnant mass to be 
small (except for large interior magnetic fields B > 10 17 G). 

A list of all simulations used to fit our model is given in 
Table |T] These results were obtained by three different groups: 
Kyoto Jl (SACRA code), UIUC JUl and the SXS collabora- 
tion J3, 8] (SpEC code). In those articles, the mass outside the 
black hole M^g 1 is measured at different times, which would 
introduce a bias in our fit. We choose to use the convention 
of Kyutoku et al. Jl, where is measured 10ms after 

merger. For this reason, the values listed in Table Udiffer from 
the masses given in the tables of J^-tll- 

Only some of the simulations listed in Table U were pub- 
lished with explicit error measurements. There is thus some 
uncertainty on the accuracy of these results. From published 
convergence tests and our own experience with such simula- 
tions, we assume that a rough estimate for the numerical errors 
AMJ,™ can be obtained by combining a 10% relative error 
and a 1% absolute error in the mass measurement, i.e. 



TABLE I: Summary of the numerical results used. When more than 
one group simulated the same set of parameters, the average value 
is used, xbh = ibh/A/bh is the dimensionless spin parameter of 
the black hole, Cns = Mns / -Rns is the compactness of the star, 
A/n°r q is the remaining mass 10ms after merger (as measured in the 
numerical simulations), and A/ NS is the baryon mass of the star. 



ID 


Af BH 


Xbh 


Cns 


NR 


Code 


Ref. 


a/ns 




1 


7 


0.90 


0.144 


0.24 


SpEC 


w 


2 


7 


0.70 


0.144 


0.05 


SpEC 


w 


3 


5 


0.50 


0.144 


0.05 


SpEC 


w 


4 


3 


0.90 


0.144 


0.35 


SpEC 


[£1 


5 


3 


0.50 


0.145 


0.15 


SpEC/SACRA 




6 


3 


0.00 


0.144 


0.04 


UIUC/SPEC 


[6,H 


7 


3 


0.75 


0.145 


0.21 


UIUC/SACRA 


£6] 


8 


5 


0.75 


0.131 


0.25 


SACRA 


[3 


9 


5 


0.75 


0.162 


0.11 


SACRA 


[3 


10 


5 


0.75 


0.172 


0.06 


SACRA 


[a 


11 


5 


0.75 


0.182 


0.02 


SACRA 


[3 


12 


4 


0.75 


0.131 


0.25 


SACRA 


[3 


13 


4 


0.75 


0.162 


0.15 


SACRA 


[3 


14 


4 


0.75 


0.172 


0.12 


SACRA 


[3 


15 


4 


0.75 


0.182 


0.07 


SACRA 


[3 


16 


4 


0.50 


0.131 


0.19 


SACRA 


[3 


17 


4 


0.50 


0.162 


0.06 


SACRA 


[3 


18 


4 


0.50 


0.172 


0.02 


SACRA 


[3 


19 


3 


0.75 


0.131 


0.24 


SACRA 


[3 


20 


3 


0.75 


0.162 


0.16 


SACRA 


[3 


21 


3 


0.75 


0.172 


0.15 


SACRA 


[3 


22 


3 


0.75 


0.182 


0.10 


SACRA 


[3 


23 


3 


0.50 


0.131 


0.19 


SACRA 


[3 


24 


3 


0.50 


0.162 


0.11 


SACRA 


[3 


25 


3 


0.50 


0.172 


0.07 


SACRA 


[3 


26 


3 


0.50 


0.182 


0.03 


SACRA 


[3 


27 


7 


0.50 


0.144 


0.00 


SpEC 


[3 


28 


3 


-0.50 


0.145 


0.01 


UIUC 


[6J 


29 


5 


0.00 


0.145 


0.01 


UIUC 




30 


4 


0.50 


0.182 


0.00 


SACRA 


[3 


31 


3 


-0.50 


0.172 


0.00 


SACRA 


[3 



is consistent with their results. We find 

a = 0.288 ±0.011 
P = 0.148 ±0.007, 



(10) 
(11) 
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for model in which tidal forces are estimated from 



(9) Newtonian physics, and 



IV. PARAMETER ESTIMATES 



a = 0.296 ±0.011 
f3 = 0.171 ±0.008 



(12) 
(13) 



A. Best-Fit parameters 

We determine the parameters a and /3 of our model (EqO 
through a least-square fit for the results of simulations 1-26 in 
TablelU Simulations 27-31, which do not lead to the formation 
of a disk, are not used directly — but we check that the model 



for the modified model M^f^ el in which the tidal forces are 
derived from the Kerr metric. 

Error estimates are easier if we rewrite the models using 
singular value decomposition (see e.g p65-75 and p793-796 
of Press et al. 1I26I1 . and references thererin), that is if we trans- 
form the basis functions of our model so that the parameters 
of the model have uncorrected errors. For example, in the 
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case of the 'Newtonian' model we have 

/i = 0.851 (3<?) 1/3 (1 - 2C NS ) - 0.525 
f 2 = 0.525 (3g) 1/3 (l-2C* NS )+ 0.851 



-Risco 
Rns 

-Risco 
-Rns 



ji.rrcm 
-'"model 



Ah + Bf 2 



The best-fit parameters A and B are then 

A = 0.323 ±0.013 
B = 0.026 ±0.001, 



(14) 



(15) 
(16) 



where the errors on A and B are independent (while the errors 
on a and (3 were strongly correlated). 



FIG. 1 : Predictions of the best-fit model (diamonds) for simulations 
1-26. The solid line represents the ideal M^ el = M^r result, 
while the error bars correspond to the estimated numerical errors 




model NS 



B. Goodness-of-fit 



C. Error Estimates 



The ability of these models to fit the numerical results 
within their errors AMj^j can be estimated through the re- 
duced x 2 



N d f 



v 26 
Z-'i— i 
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model 
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where N * = 26 — A^ params = 24 is the number of degrees 
of freedom, and the index i refers to the ID of the numerical 
simulations. The 'Newtonian' model M™^ el and the 'Kerr' 
model M, 



model 



are equally good fit to the data, with \ = 
0.98 and \ 2 = 0.96 respectively. By comparison, the best-fit 
results for model M modcl (in which we do not correct dtidai by 
the factor [1 — 2Cns]) has a much larger \ 2 = 4.04. Adding 
a constant term 7 to either M^™^ or M™odei leads to x 2 = 
1.00. 

A comparison between the simple model M£f™[ el and the 
numerical results is shown in Fig.Q] in which we plot AJ^r 11 
as a function of MJ^ el for simulations 1-26. We can see that 
the difference between the modelled and measured masses is 
generally smaller than the errors expected from Eq. [9] The 
main exception is the large remnant mass observed in case 
4. We suspect that our model, which assumes that the rem- 
nant mass scales linearly with -Risco an d ^tidab breaks down 
for remnant masses greater than about 20 — 25% of the neu- 
tron star mass. A non-linear relation between these distances 
and the remnant mass might perform better in that regime, but 
more numerical simulations are required to test that hypothe- 
sis. 

The more complex model A/' 1 ™ lcl offers very similar re- 
sults: for cases 1-26, the worst disagreement between the 
models is 0.008M^ S (for case 3) while their rms difference 
is 0.004M& S . 



Estimating the error in the mass predictions of our model 
from the statistical errors in the parameters a and j3 is likely 
to be misleading. Differences between the numerical results 
M^™ 1 and the predictions of the model AT™ °j el come from 
multiple sources: the numerical error AA/^™ of course, but 
also a physical spread of the exact mass remnants around the 
predictions of the model. A part of that spread at least should 
be due to differences in the outcome of BHNS mergers for 
binaries with the same parameters (A/bhi-^ns>Xbh>Cns)> but 
different equations of state (i.e. neutron stars with the same 
radius but a different internal structure). But more generally, 
it is unlikely that the simple equations used here can perfectly 
represent the complex dynamics of a BHNS merger. 

From the fact that we measured \ 2 ~ 1, we know that the 



errors Af™ c| 



MJj™ are compatible with a gaussian distri- 



bution of variance AA/^™ 1 - This is already indicative of the 
likely existence of a non-zero physical spread around the re- 
sults of the model. The estimated numerical errors AA/^™ 
are indeed more of an upper-bound on the errors in the simu- 
lations than the width of an expected gaussian distribution. In 
the absence of a difference between the real physical outcome 
of a merger and the output of the model, we would thus ex- 
pect x 2 t° be lower than 1 . How much of the measured errors 
Afmodci — MJ,™ comes from numerical errors and how much 
from actual differences between the model and the physical 
reality is hard to determine, especially considering that the nu- 
merical errors are not well known. A more cautious approach 
to estimate the uncertainty in the model is thus to consider 
AMJ,™ as a conservative upper bound on the variance of the 
error in the model itself. 

Fig. shows contours of M^ el = 0.1M£ S for various 
neutron star compactness. The general features of this plot 
are not surprising: the formation of massive disks is known to 
be favored by low mass ratios, high black hole spins and large 
neutron stars. But our model allows for the determination of 
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FIG. 2: M£J™i 8l = O.IA/Js contours for, from top to bottom, neu- 
tron star compactness Cns = 0.22,0.18,0.155,0.135 (i.e. R^s re 
9.5, 11.5, 13.5, 15.5km for Mns = 1.4M©). For each compact- 
ness, we have M™™± el > 0.1Mn S above the plotted contour. The 
shaded regions encompass the portions of phase space for which 
M£™dei = O.lMSs ± &Mnr- SGRBs are extremely unlikely to 
occur below the green region (Cns = 0.155). Note that the scale is 
chosen in order to zoom on the high-spin region. 




the region of parameter space in which a certain amount of 
matter will remain available at late times with fairly high ac- 
curacy: at least within the spread AMJfg re 0.02M© or, if 
we consider a measurement of M Tem as a way to determine 
the radius of a neutron star, within Ai?NS ^5 0.5km. 

Another important issue is the validity of the model outside 
the parameter range currently covered by numerical relativ- 
ity. It is indeed possible that larger errors will be found for 
more compact neutron stars (Cns > 0.18) or larger mass 
ratios (A/bh > 7A1ns)- However, given that our model 
fits the numerical data over a fairly wide range of parame- 
ters, and is derived from the physics of tidal disruption, it 
is likely to give decent results over most of the astrophys- 
ically relevant parameter space — with the notable excep- 
tion of configurations leading to very large remnant masses 
M rcm > o.20-0.25M NS (i.e. for nearly-extremal black hole 
spins and low mass ratios), and probably of the asymptotic 
regime xbh — > 1 where scalings valid in the range \ = 0—0.9 
might break down. 



V. RESULTS 

The models described in the previous sections can be used 
to easily approximate the region of parameter space in which 
disruption occurs, or in which a certain amount of mass will 
remain available at late times. Such predictions are particu- 
larly important when trying to determine which BHNS merg- 
ers are likely to lead to short gamma-ray bursts (SGRBs). 
BHNS mergers ending with the formation of a massive ac- 
cretion disk might power SGRBs. In the absence of a disk, 
on the other hand, the only observational signatures of BHNS 
mergers are their gravitational wave emissions, as well as po- 
tential electromagnetic or neutrino precursors (see e.g. Tsang 



et al. il). 

The minimum remnant mass required to get SGRBs is cur- 
rently unknown, and is likely to vary across the parameter 
space: the fraction of the remnant mass which, at any given 
time, is in a long-lived accretion disk around the black hole 
(as opposed to the tidal tail or unbound ejecta) is by no means 
a constant, nor are the physical characteristics of that disk. 
We know, for example, that at high mass ratios a larger frac- 
tion of the mass is initially in an extended tidal tail than for 
lower mass black holes Q. Nevertheless, M^™[ el is already 
a useful prediction, providing a good estimate of the amount 
of material available for post-merger evolution. Additionally, 
any configuration for which M^™[ eI = can be immediately 
rejected as a potential SGRB progenitor. 

Predictions for the mass of neutron star material remaining 
outside the black hole at late times are detailed in Figs. [3]|5] 
in which we plot contours of the remnant mass as a function 
of the mass ratio and black hole spin. Each figure correspond 
to a different neutron star compactness, covering the range 
of radii expected from the theoretical results of Hebeler et 
al. l22ll . Experimental measurements of neutron star radii are 
still fairly difficult, but studies of bursting X-ray binaries by 
Ozel et al. [28-30] tend to favor the lower range of potential 
radii (i?NS re 9 — 12 km). Steiner et al. ll23ll . after reassessing 
the errors in the measurement of neutron star radii from X-ray 
bursts, derived a parametrized equation of state which takes 
into account both the astrophysical measurements and results 
from nuclear theory. They predict that .Rns ~ 11 — 12 km 
for Mns = IAMq. We can thus consider Fig.|3]and Fig.|5]as 
bounding the range of potential neutron star radii, while Fig. [4] 
is around the most likely neutron star size. 

The strong dependence of the remnant mass in the radius 
of the star is particularly noteworthy. In the most likely astro- 
physical range of mass ratios (q ~ 5 — 10), remnant masses 
M rem = O.IMns can be achieved for moderate black hole 
spins xbh ~ 0.7 — 0.9 if we consider neutron stars with 
C NS = 0.155 (R NS re 13.5km), as in Fig. [3] But at the 
other end of the range of potential neutron star radii, for 
Cns = 0.22 (i?NS ~ 9.5km), the much more restrictive 
condition xbh ~ 0.9 — 0.999 applies(Fig. [5). For a neu- 
tron star in the range of compactness favored by Steiner et 
al. JH (C NS = 0.18, or R NS re 11.5km), keeping 10% 
of the neutron star material outside the black hole requires 
spins xbh ~ 0.8 — 0.97, an already fairly restrictive condi- 
tion (Fig. ID . 

This naturally implies that the rate of SGRBs produced as 
a result of BHNS mergers is very sensitive to the equation of 
state of nuclear matter, and in particular to the size of neutron 
stars. Determining that rate is unfortunately impossible with- 
out knowledge of the number of BHNS mergers, and of the 
distributions of black hole spins and mass ratios. Additionally, 
M lcm is not a fully reliable predictor of the ability of a given 
BHNS binary to power a SGRB. Nonetheless, the importance 
of the equation of state can be fairly easily understood by 
simply computing the area of the region above a certain con- 
tour of MJf™Lj for various values of Cns- Let us define 
X C (M, Cns> q) as the critical spin above which M™^ ol > M 



6 



'model = (0,0.05, 0.1, 0.15, 0.2)Af£ s for a 
star of compactness Cns = 0.155 (i?NS ~ 13.5km for Mns = 
IAMq). The shaded regions correspond to portions of parameter 
space in which no matter remains around the black hole (bottom/red), 
more than 0.2Af£s remains and massive disks should be the norm 
(top/green), and an intermediate region in which lower mass disks 
will form (center/blue). Note that the scale is chosen in order to 
zoom on the high-spin region. 
0.95 1 



H 0.8 




FIG. 4: Same as Fig.|3] but for Cns = 0.18 (i? N s ~ 11.5km). 
0.99 r 



0.95 



0.9 




FIG. 5: Same as Fig.0 but for Cns = 0.22 (7? N s « 9.5km). 
0.999 r 



0.995 - \M=0.2M 



0.99 - 




0.95 



TABLE II: Fraction 4>(M, Cns) of the parameter space within q — 
5 — 10, xbh = — 1 for which A/modci > M for various neutron 
star compactness Cns and critical masses M. 



Cns 


<^(0,Cns) 


<^(0.1M N s,Cns) 


<^(0.2M ns ,Cns) 


0.135 


0.46 


0.30 


0.16 


0.155 


0.29 


0.17 


0.07 


0.180 


0.16 


0.08 


0.02 


0.220 


0.05 


0.01 


0.00 



and 



0(M,C NS ) = 



J 5 10 [l 



X c (M,C NS ,q)}dq 



(18) 



Then, </>(M, Cns) represents the fraction of binaries with 
mass remnants greater than M assuming that the distributions 
of mass ratios and spins are uniform within the q = 5 — 10 
and xbh = — 1 range respectively. As we decrease the 
size of the neutron star from Cns = 0.155 to Cns = 0.22, 
Table [TT] shows that we go from about 20% of the parame- 
ter space in which significant disks are possible to about l%! 
This does not necessarily mean that SGRBs are impossible for 
Cns ~ 0-22 — but certainly indicate that they would occur 
in a non-negligible fraction of BHNS mergers only if quasi- 
extremal spins are the norm. 

These predictions can be compared with previous estimates 
for the disk mass obtained by Pannarale et al. lfl9ll . In lfl9tl . 
an estimate of the mass available at late times is obtained 
through the use of a toy-model studying the tidal forces act- 
ing on the neutron star, which is represented by a tri-axial el- 
lipsoid. If the qualitative features of both models are fairly 
similar, we find that the use of more recent numerical results 
(and, in particular, the availability of mass estimates for larger 
mass ratios) causes significant differences in their predic- 
tions. Specifically, the results presented here indicate that it 
is a lot more difficult to create massive accretion disks at high 
mass ratios than what that previous model indicated. On the 
other hand, at low mass ratios q ~ 3 our model predicts signif- 
icantly higher final masses — probably because many of the 
numerical simulations used to calibrate lfl9tl are now known 
to have underestimated the amount of matter remaining out- 
side the black hole. As opposed to ifl^l . our model is unlikely 
to capture the behavior of BHNS mergers with q ~ 1, when 
finite-size effects begin to make it more difficult to form mas- 
sive disks. 

We can also revisit the condition derived by Taniguchi et 
al. II18[| for the parameters allowing disk formation in the case 
of non-spinning black holes. Requiring M™£^ ei > for 
Xbh = is equivalent to imposing an upper bound on the 
neutron star compactness, 



°NS <~ 



1 



2 + 2.14g 2 / 3 ' 



(19) 



For massive black holes (q > 7), this result is very sim- 
ilar to what Taniguchi et al. 1180 found by studying quasi- 
equilibrium sequences of BHNS binaries. At low mass ratio, 
our results are less favorable to tidal disruption and disk for- 
mation. This is in agreement with the numerical simulations 
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7M NS (M BH w 10M o ). 
(0,0.05, 0.1, 0.15, 0.2)M& S . The shaded regions are as in Fig. [3] 
and the neutron star radius (top scale) is computed assuming a star of 
ADM mass M NS = 1.4M . 



R ns 0™) io 



12 



13 



14 



1 1 

Massive Disk / SGRB 




performed by Kyutoku et al. U, even though the results for 
low mass ratio, non-spinning BHNS mergers published in 121 
were not taken into account when fitting our model. 

Current population synthesis models estimate the peak of 
the distribution of black hole masses in BHNS systems to be 
around M BH ~ 10M Q , or M BH ~ 7M NS HQ. Fig. gof- 
fers clearer information on the behavior of BHNS systems in 
that regime. We see that no disk can form for xbh < 0.9 un- 
less i?NS > 10.5km. That condition becomes i?NS > 12km 
if we require at least 0.1M NS outside the black hole at late 
times. Results for BHNS binaries with higher black hole spins 
(Xbh — > 1) should of course be considered with caution: in- 
deed, no mergers of BHNS binaries have been published for 
Xbh > 0.9 or Cns > 0.18, and such simulations would 
be required to rigorously test the accuracy of these predic- 
tions in extreme regions of the parameter space. Nonethe- 
less, our model indicates that quasi-extremal spins are at least 
a necessary condition for the formation of massive disks for 
M B h ~ 10M Q and R NS < 12km. 

The minimum spin requirement for massive disk forma- 
tion across the parameter space of BHNS binaries is shown 
in Fig. [7j in which the black hole spin needed to keep 10% 
of the neutron star mass outside the black hole at late times 
is plotted. Figs. [6] and [7] both indicate the existence of an 
extended region of parameter space (Cns ~ 0.18 — 0.22, 
Xbh ~ 0.9 — 1) which is likely to be astrophysically rele- 
vant but remains numerically unexplored, and in which the 
outcome of BHNS mergers varies significantly. 



VI. CONCLUSIONS 

We constructed a simple model predicting the amount of 
matter remaining outside the black hole about 10ms after a 
BHNS merger, based on comparisons between the binary sep- 



J modoi — O.IMns for varying black hole spins 
Xbh = (0.5,0.7,0.9,0.99). The grey region contains spins above 
the maximum value reached by numerical simulations. 




aration at which the neutron star is expected to be disrupted 
by tidal forces from the black hole and the radius of the inner- 
most stable circular orbit around the hole. We show that the 
model can reproduce the results of recent general relativistic 
simulations of non-precessing, low-eccentricity BHNS merg- 
ers within a few percents of the total mass of the neutron star. 
The simplest best-fit model is 



JW NS 



0.288 3 



M 



BH 



Mi 



NS 



1/3 



1 - 2- 



.M 



NS 



V 



Rns J 



0.148 



-Risco 
-Rns 



(in units in which G = c = 1). 

These mass predictions should be valid at the very least 
within the range of parameters currently covered by numer- 
ical simulations (Mbh = 3 — 7Mns ; -Rns = 11 — 16 km, 
qbh / Mbh = — 0.9), and are likely to remain fairly accurate 
within most of the astrophysically relevant parameter space. 

Using this model, it becomes easy to estimate the region 
of parameter space in which large amounts of matter remain 
outside the black hole for long periods of time. This is of par- 
ticular importance when studying whether BHNS mergers can 
result in short gamma-ray bursts. Our results show the strong 
dependence of the remnant mass in the radius of the neutron 
star: whether the equation of state of neutron stars is at the soft 
or stiff end of its potential range could easily translate into an 
order of magnitude difference in the rate of gamma-ray bursts 
originating from BHNS mergers. It is also quite clear that 
high black hole spins are likely to be a prerequisite for the 
formation of massive disks. Neutron stars in the middle of the 
theoretically allowed range of radii (-Rns ~ 11.5km) require 
spins obh /Mbh ^ 0.8 for about 10% of the neutron star ma- 
terial to remain outside the hole, while quasi-extremal spins 
are necessary for the most compact stars. 

The validity of our model is currently limited to black hole 
spins aligned with the orbital angular momentum and rem- 
nants below ~ 20 — 25% of the neutron star mass, due to the 
lack of numerical data available for precessing binaries and 
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high mass remnants. Extending the model to cover these inter- 
esting parts of the parameter space would certainly be useful, 
but would require a large number of computationally intensive 
simulations to be performed (particularly to cover misaligned 
black hole spins). A few additional simulations using high 
mass ratios or small neutron star radii together with relatively 
large spins (xbh ^ 0.9) would also be extremely helpful, al- 
lowing better estimates of the errors in the model for binary 
parameters which are astrophysically relevant but have never 
been considered by numerical relativists. 

The extreme simplicity of these models should make them 
useful tools to obtain cheap but reliable estimates of the results 
of BHNS mergers across most of the astrophysically relevant 
parameter space, as well as to help determining which numer- 



ical simulations to perform in order to study given physical 
effects. 
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